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Abstract 

By proposing a statistic for the detection of secondary (Doppler) peaks in the CMBR 
power spectrum, the significance level at which such peaks may be detected are computed 
for a large range of model CMBR experiments. In particular, we investigate experimental 
design features required to distinguish between competing cosmological theories, such as 
cosmic strings and inflation, by establishing whether or not secondary peaks are present 
in the CMBR power spectrum. 

1 Introduction 

Experimental measurement of the Doppler peaks' positions and heights would fix at least some 
combinations of cosmological parameters (e.g. Hq, Qq etc.) which are left free in inflationary 
models Q. Furthermore, as shown by [@, HJ, the absence of secondary Doppler peaks is a robust 
prediction for cosmic strings, although this may or may not be the case for textures [B, [|]. 
Therefore it appears that even determining whether or not there are secondary Doppler peaks 
offers an important alternative motivation for measuring the CMBR power spectrum. 

We address this issue by proposing a statistic for detecting secondary oscillations, and 
studying how it performs for various models, using different experimental strategies. The 
results are encoded in a detection function E which indicates to within how many sigmas we 
can claim a detection of secondary oscillations, given a particular model and experiment. 

We apply the statistic to both the standard CDM scenario and an open CDM model which 
is tuned to confuse inflation and cosmic strings in all but the existence or otherwise of secondary 
oscillations. For a wide range of experiments we allow the beam size, sky coverage, and detector 
noise to vary, and use this framework to compute the detection function for secondary peaks. 



2 Power spectrum estimation from real observations 



There are several factors affecting how well one can measure the CMBR power spectrum from 
real observations, which we now discuss. 

(i) Distortion of the underlying spectrum due to the finite size of the observed field. For 
(square) fields of size L > 4 degrees (suitably windowed with a cosine bell or Hann window) 
this is not a severe problem for detecting secondary peak structure. 

(ii) Cosmic/sample variance, which places constraints on the minimum sky-coverage neces- 
sary to achieve a given accuracy. Roughly speaking, if / s is the fraction of sky observed, then 
* 2 (C e )/Ci*l/(£f s ). 

(iii) Instrumental noise, which we shall assume is uncorrelated for simplicity, and charac- 
terised by o"pi x , the rms pixel noise, and fi P i x , the area of a pixel. If we consider the most general 
case where only a fraction / s of the sky is mapped, then for a detector of fixed sensitivity, and for 
a fixed total observing time, then by varying f2 pix and / s the quantity w^ 1 = crp ix f2 pix (47r / / s ) , 
remains constant, and is therefore an important qualifier for noise on maps obtained using 
different scanning strategies. 

(iv) Diffuse foreground emission, which can severely hamper the measurement of CMBR 
anisotropies. A discussion of these foreground components, and the regions of frequency/multipole 
space in which each dominates, is given by ||. The main components of this foreground are 
Galactic dust, synchrotron and free- free emission. Algorithms for separating these components 
from the CMBR signal are discussed by || and [0. Typically the errors associated with the 
separation process are of a similar magnitude to the average errors on an individual frequency 
channel due to instrumental noise alone, but details depend on the separation algorithm used. 

(v) Point sources, which cannot be removed from spectral information alone. This requires 
the identification of the sources by higher-resolution observations at a frequency close to that 
of the CMBR observations, with sufficient flux sensitivity to indentify all point sources down 
to some flux limit roughly equal to the instrumental noise of the CMBR observations. We 
note here that although it is generally believed that point source contamination becomes less 
important as the observing frequency increases above about 100 GHz, there is no direct evidence 
for this. Moreover, even the population of radio point sources at frequencies above about 10 
GHz is rather uncertain, and it may be inadvisable to rely on low frequency surveys such as 
the 1.5 GHz VLA FIRST survey || to subtract point sources from CMBR maps made at much 
higher frequencies. 

3 Observing Doppler peaks in standard CDM 

The idea is to apply to a particular model a statistic sensitive only to the existence or absence 
of secondary oscillations in the power spectrum. In this section we consider the power spectrum 
predicted by the standard inflation/CDM scenario with Q = 1, h = 0.5 and Oj, = 0.05 (which 
we shall call sCDM). 

To this end we first compute the average broad band power Cj in each of several equally 
spaced-bins, denoted by horizontal bars in Fig. [l| We then infer the convexity Ci of the spectrum 
at each bin position (apart from the first and last bins) from d = (Cj_i + Cj + i)/2 — C{. These 
convexities are all negative if there are no secondary peaks, but alternate in sign for sCDM. If 
the overall error in Ci is a 2 (Cj) then one can define an oscillation detection function as 
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Figure 1: The ensemble average power spectrum predicted by sCDM (solid line), and the en- 
semble average power spectrum for a square observing field of size 10 degrees after windowing 
with a cosine bell (dashed line). The points indicate the corresponding power spectrum esti- 
mates for a particular realisation in each of the bins indicated by the horizontal error bars. The 
vertical error bars indicate the theoretical cosmic/sample variance in the absence of noise. 



for i = 2 and i = 4, which tells us to within how many sigmas we can claim a detection 
of secondary peaks. The method for computing the estimates Cj of the power spectrum in 
each bin, and their associated standard errors, taking into account limited sky-coverage and 
instrumental noise, are discussed in detail in || and [10|. 

From Fig. |1] we see that the first dip in the sCDM power spectrum is more easily detected 
than than the second one, a situation only exacerbated by finite resolution and the presence of 
instrumental noise. Therefore we shall confine ourselves to considering the detection function 
E 2 , which from now on we refer to simply as E. 

The detection function E = E(L, 6b, w -1 ), where L is the linear size of the observed (square) 
field [the all-sky limit can be recovered by setting L 2 = 4it sr « (202 deg) 2 ], 6b is FWHM of the 
observing beam, and w~ l is the noise level discussed above. This function is plotted in Fig. |2| 
for the low noise case w^ 1 = (25/iK) 2 (deg) 2 and the high-noise case w^ 1 = (60/iK) 2 (deg) 2 . 

For any beamsize there is a maximum sky coverage beyond which the detection is not 
improved. If anything the level of the detection decreases, but typically not by much. The 
ideal scanning strategy is then defined by a line Li(6b) which intersects the contours of E at the 
lowest L-value at which a plateau has been achieved in the detection function. The significance 
of the detection obtained for an ideally scanned experiment depends on the beam size. For 
example, in the low-noise case, if 6b = 0.6°, the ideal coverage is a patch of Lj(0.6°) = 5 degrees, 
which results in a 3-sigma detection. If 6b = 0.5°, on the other hand, an 8-sigma detection can 
be obtained with Li = 35 degrees. The detection provided by an optimally scanned experiment 
increases at first very quickly as the beam is reduced below 6b = 0.6° (from 3-sigma at 6 = 0.6° 
to 33-sigma at 6b = 0.2°). By reducing 6b from 0.2° to zero, however, the detection is only 
increased by 2-sigma (from 33 to 35). For this level of noise the maximal detection is 35 sigma 
and is achieved with 6b < 3' and all-sky coverage. For low noise levels all-sky coverage is never 
harmful, but it is the beamsize that determines how good a detection can be achieved, and how 
much sky coverage is actually required for an optimum level of detection. 
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Figure 2: (Left) Low noise u> _1 = (25fiK) 2 (deg) 2 contours of the £ function. Sky coverage L 
varies between 5 degrees and all sky (L = 202 degrees), and the beamsize FWHM between 
and 1 degrees. (Right) High noise w -1 = (60[iK) 2 (deg) 2 contours of the £ function. 

For noise levels of the order w" 1 = (25fiK) 2 (deg) 2 the overall picture is always as in Fig. ^|. 
In particular, there is always a top contour indicating the maximal detection allowed by the 
given noise level. The maximum £ is always achieved with infinite resolution, but one falls 
short of this maximum by only a couple of sigmas if 6^ m 0.1°. If the noise is much smaller than 
this, however, the summit of £ is beyond L = 202°. For w^ 1 = (15[iK) 2 (deg) 2 , for instance, 
all-sky coverage becomes ideal for any dt, < 0.3°. 

If, on the other hand, the noise is much larger than w~ l = (25fiK) 2 (deg) 2 then the £ 
contours are qualitatively different, as shown in Fig. |2|for w~ l = (60fiK) 2 (deg) 2 . The beamsize 
is now a crucial factor. A beamsize of 6^ = 0.5° would provide a 3-sigma detection (with 
Li = 10°), but reducing the beamsize to about 6^ = 0.4° improves the detection to 6-sigma 
(with Li = 20°). It is also clear from the figure that, for high noise levels, forcing all-sky 
coverage dramatically decreases the detection. 

4 Open CDM models and cosmic strings 

We may repeat the above analysis for different cosmological models. We therefore consider the 
case of maximal confusion between inflation/CDM and cosmic string scenarios by comparing a 
cosmic strings model with a CDM model for which the main peak in the power spectrum has 
the same position and shape (but the latter exhibits secondary peaks) . For definiteness we have 
chosen a CDM theory with a flat primordial spectrum, = 0.3, ho = 0.6, and flbh^ = 0.02. 
We shall call this theory stCDM, the CDM competitor of cosmic strings. As before we simply 
study the first dip detection function of stCDM, and then take this detection function as a 
cosmic string rejection function. 

In Fig. ^ we show the angular power spectrum of stCDM (solid line) and a possible power 
spectrum for cosmic strings (dotted line). We then simply repeat the same exercise as in the 
previous section to obtain the detection function of the first dip of stCDM. The results are 
shown in Fig. |] for the same noise levels as before. 

Overall we see that in signal dominated regions the detection is much better for stCDM than 
for sCDM. This is because features at higher I have a smaller cosmic/sample variance (which is 
proportional to !/£)■ It can be checked that the cosmic/sample variance limit, obtained with a 
single-dish experiment with no noise, is now E w 197L 2 /(Air) (as opposed to E ~ 77L 2 /(47r) for 
sCDM). Even in the presence of noise, wherever the signal dominates, the detection is better 





10 



100 

I 



1000 



Figure 3: The angular power spectrum of stCDM (solid line) and one possible cosmic string 
scenario (dotted line). The dashed line is the ensemble average stCDM power spectrum as 
sampled by an experiment with a field size L = 2. The points indicate the average power in 
each bin for stCDM. The horizontal errorbars denote the width of the bins, and the vertical 
errorbars show the sample variance of the power estimates for such an experiment assuming no 
instrumental noise. 

for stCDM. However, in noise-dominated regions the behaviour of the detection function for 
stCDM and CDM is very different. 

The signal-dominated region is greatly reduced in stCDM. Much smaller beamsizes 9b are 
now required for any meaningful detection. As shown in Fig. |], one would now need 6b < 0.3° 
and 6 b < 0.25°, for noises w^ 1 = (25/iif) 2 (deg) 2 and w^ 1 = (60/ii^) 2 (deg) 2 respectively, in order 
to obtain a reasonable detection. Again one can plot an ideal scanning line in the beam/coverage 
sections defined by a fixed noise u> -1 . The ideal sky coverage is much smaller for stCDM than 
for sCDM. In general the contours of £ for stCDM compared to sCDM are squashed to lower 
6b, lower L, and achieve higher significance levels, with steeper slopes. Following an ideal 
scanning line for any fixed w" 1 one reaches a maximal detection allowed by the given level of 
noise, which is always better for stCDM than for sCDM. This maximal detection is normally 
obtained with a small sky coverage, and infinite resolution. Nevertheless, one falls short of this 
maximum by only a few sigma if the resolution is about 9b = 0.05° — 0.1°. From Fig. |], for 
w^ 1 = (25fiK) 2 (deg) 2 , one may now obtain a maximal 43-sigma detection for an ideal scanning 
area of L = 65 degrees. If 9b = 0.1° a 36-sigma detection is still obtained. We also see that 
a beamsize of 9b < 0.25° is required to obtain a 3-sigma detection (with L = 4 degrees), and 
a 10-sigma detection can be achieved only with 9b ~ 0.15° (and Li = 18 degrees). All-sky 
coverage for an experiment targeting stCDM is generally inadvisable, and it would only be 
optimal for the extremely low level of noise w~ l < (11 fiK) 2 (deg) 2 . 

5 Conclusions 

The results obtained here are useful for future CMBR projects in two different ways. Firstly 
they allow the choice of an ideal scanning strategy (choice of resolution and sky coverage) for 
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Figure 4: (Left) Low noise w~ l = (25/ifT) 2 (deg) 2 contours of the detection function X for 
stCDM. (Right) High noise w~ l = (60/iK) 2 (deg) 2 contours of the detection function. 

detecting secondary Doppler peaks, given observational constraints such as the instrumental 
noise level and the total observing time. Secondly, one may compute the expected value of 
the detection, assuming ideal scanning, as a function of these parameters. This provides lower 
bounds on experimental conditions for a meaningful detection as well as an estimate of how 
fast detections will improve thereafter. 

These results also indicate that in order to study Doppler peak features for sCDM, depending 
on the noise levels, a large sky coverage might be desirable, even for a resolution of about 
6b = 0.4° — 0.5°. If, however, one is instead to test the high-/ opposition between low Q CDM 
and cosmic strings, then a rather higher resolution is required. Furthermore, in this context, 
all-sky scanning is not only unnecessary, but in fact undesirable. 
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